Many-body models for molecular nanomagnets 
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We present a flexible and eflFective ab-initio scheme to build many-body models for molecular 
nanomagnets, and to calculate magnetic exchange couplings and zero-field splittings. It is based on 
using localized Foster-Boys orbitals as one-electron basis. We apply this scheme to three paradig- 
matic systems, the antiferromagnetic rings Crs and CryNi and the single molecule magnet Fe4. In all 
cases we identify the essential magnetic interactions and find excellent agreement with experiments. 
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Clusters made of a finite number of interacting spins 
are ideal test beds to investigate fundamental issues in 
quantum mechanics. One of the first physical realizations 
are molecular nanomagnets (MNMs) , molecules contain- 
ing a core of d or / ions, whose spins are coupled by 
magnetic interactions; MNMs form crystals which behave 
like an ensemble of identical and almost non-interacting 
magnetic units. During the last years sophisticated ex- 
periments and targeted research activities have unveiled 
a variety of fundamental quantum phenomena and poten- 
tial technological applications of MNMs [ll-Q- The two 
most promising classes have been identified in the Single- 
Molecule Magnets (SMM), like Mni2 Q, Fesfl and Fe4 
P, SjS, and the Antiferromagnetic Rings (AFR), like 
Crg lia [H, shown in Fig. [H CrrNi H, [H [ll and Fee 
While SMM have opened the perspective of storing 
information in single molecules and building high-density 
magnetic memories, AFR are of great interest in the field 
of quantum information processing 



At the synthetic level, thanks to the huge progresses 
made in the last years in coordination chemistry, it is now 
possible to reach a high degree of control on the molecular 
structure and on the topology of magnetic interactions. 
At the theoretical level, one of the main obstacles to fur- 
ther progress remains the lack of a flexible and systematic 
approach to build ab-initio system-specific models for the 
magnetic interactions; such models should describe on 
the same footing chemistry and many-body effects within 
the partially filled d or / shells of the magnetic ions. 

MNMs are typically described through Heisenberg-like 
spin Hamiltonians. If the form of the spin Hamilto- 
nian is known, the magnetic couplings can be in princi- 
ple extracted from total-energy density-functional theory 
(DFT) calculations for different spin configurations 16- 
|l9[. This approach can become unpractical if many pa- 
rameters have to be determined, as e.g. in heterometallic 
compounds or anisotropic SMM; furthermore, subtle in- 
teractions, which could greatly influence, e.g., the relax- 
ation dynamics, can be easily overlooked. An alternative 
consists in computing the couplings via energy variations 
at small spin rotations (20| . However, at a more funda- 
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FIG. 1: The Crg antiferromagnetic ring and its Cr d crystal- 
field orbitals for site 1, in order a) — > e) of increasing energy. 
Red (blue) orbital lobes are positive (negative). The covalent 
p tails on the neighboring ligands carry the information on 
the molecular structure and are crucial for the magnetic ex- 
change couplings. We define z as the axis perpendicular to the 
ring and pointing outwards. All Cr sites are approximatively 
equivalent. H atoms are not shown for clarity. 



mental level, a common problem of all these approaches 
is that the most used DFT functionals (the local-density 
approximation (LDA) and its simple extensions), do not 
properly describe strong correlation effects in open d- 
or /-shells, while LDA-\-U or hybrid functionals include 
them only at the static mean-field level. Recently, it has 
been suggested that Hubbard-like models could be more 
appropriate 2lJ, |22| , but an efficient and flexible scheme 
to calculate the parameters of such models ab-initio, no 
matter the complexity of the system, has not been im- 
plemented so far. 

In this Letter we show that this can be achieved by 
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FIG. 2: Structure of Fe4 and Fe d crystal-field orbitals local- 
ized on the central Fe (-D3 symmetry): The ai singlet (a), and 
the two e doublets (&, c) and (d, e), in order of increasing en- 
ergy. Heisenberg couplings: F^'^ (full line) and F^''^ (dashed 
line). We define z as the axis perpendicular to the triangle, 
and pointing outwards. H atoms are not shown for clarity. 



using localized Foster-Boys orbitals [23| as one-electron 
basis to construct molecule-specific generalized Hubbard 
models. We use the constrained local-density approx- 
imation (cLDA) scheme [l^] to calculate the screened 
Coulomb interactions in such a Foster-Boys basis. We 
obtain the spin Hamiltonians systematically by using a 
canonical transformation 25| to eliminate charge fluctua- 
tions, without any a-priori assumption on the form of the 
final spin Hamiltonian. We implement this scheme in the 
NWChem quantum-chemistry code [2^, and apply it to 
three prototype molecules, representative of the two main 
classes of MNMs: the AFRs Crg (Fig. [U and CryNi and 
the SMM Fe4 (Fig. These systems have been exten- 
sively investigated experimentally [sl-fiol [l2. 27 1 and the 
magnetic exchange couplings are now well known. In all 
cases, we find excellent agreement with experiments, and 
identify the microscopic mechanisms that lead to the em- 
pirical spin models commonly adopted to describe them. 

The procedure we adopt is the following. First we per- 
form LDA calculations for the experimental structures 
reported in Refs. in this step we use as basis a 

triple-zeta valence set of gaussians. Next, we identify the 
transition-metal o?-like molecular orbitals; by means of 
Foster-Boys localization [H] , we construct a set of local- 
ized orbitals, which are centered on the transition-metal 
ions and span such c?-like states; using these orbitals as 
basis we build the corresponding generalized Hubbard 
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Here c|„jg. (Cj^^.) creates (annihilates) an electron with 
spin a in the Boys orbital m at site i. The parameters 
^mm' ^^'^ hopping integrals (i ^ i') or the crystal- 
field matrix (i — i'), while 11^^^,^, are the screened 
Coulomb integrals (Tab. |l|. The term Hbc is the dou- 
ble counting correction, which removes the part of the 
Coulomb interaction already included and well accounted 
for in the LDA; A,; is the spin-orbit coupling. The re- 
sults presented in this work are obtained using the ro- 
tational invariant form of the Coulomb vertex, includ- 
ing spin-flip and pair hopping terms but (for simplicity) 
no Coulomb anisotropy; thus all Coulomb parameters 
can be expressed as a function of the averaged screened 
Coulomb couplings t/*'' and J*'' We determine the 
latter by using the cLDA [i^l approach in the Foster-Boys 
basis, keeping the basis frozen in the self-consistency 
loop. For Hdc we adopt the common (30| expression 
^^DC = ^E.t^'''"^K-l)-3E.'/''X(H-l)'^here 
is the number of d electrons at site i. For homonu- 
clear systems Hbc amounts to a shift of the d levels 
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Fe4 


CryNi 




5.98 


5.22 


6.32 (Ni) 


[/2.2 


5.98 


5.03 


5.98 (Cr) 




0.26 


0.24 


0.21 (Ni) 


j2,2 


0.26 


0.22 


0.26 (Cr) 


Al 


16.5 


34.3 


33.5 (Ni) 


A2 


16.5 


37.0 


16.5 (Cr) 



TABLE I: Top: Crystal-field energy levels and hopping inte- 
grals for Crg. The latter are given in the basis of crystal-field 
orbitals (in eV), and for sites i = 1 and i' = 2 in Fig. [1] The 
energy of the Fermi level is set to zero. Bottom: Screened 
Coulomb integrals and J*'' obtained via cLDA. Sites i 
are defined in Fig. \T\ and Fig. [5] 



3 





I 




10 12 

U (eV) 



FIG. 3: Crg: Calculated super-exchange coupling (in meV) 
between nearest neighbors [sj], Tg^"^^, as a function of the 
Coulomb parameters J and U . Calculations are performed 
using the rotationally-invariant Coulomb interaction. The 
cLDA values of the screened Coulomb integrals are f/ ~ 6 eV 
and J ~ 0.26 eV. Around these values Eg'^^^ > (AFM). 



and can be incorporated in the chemical potential; in 
the case of CrrNi, instead, the shift due to Hj^c has to 
be taken into account explicitly. Finally, we extract the 
spin-orbit coupling by comparing the one-electron part 
of Hamiltonian ([1]) obtained with and without spin-orbit 
interaction. Once we have obtained the parameters of 
the Hubbard model (Tab.|I|), by using a canonical trans- 
formation, we eliminate charge fluctuations and derive 
the corresponding low-energy spin model. In this step, it 
is convenient to work in the basis of crystal-field states, 
obtained by diagonalizing the on-site matrix i^'^/; we 
denote their energies £„ with e„ < En+i- The crystal- 
field states are shown in Fig.[T]for Crg and Fig.[2]for Fe4. 
At all sites but the central Fe in Fe4, the environment 
of the magnetic ion is approximatively octahedral; thus 
the crystal-field orbitals split into a lower energy t2g-like 
quasi-triplet and a 1 — 2 eV higher energy e^-like quasi- 
doublet. The central Fe site of Fe4 has D3 symmetry; 
its crystal-field levels (Fig.[2|) split into a ai ground state 
and two excited e doublets, ^ 0.6 eV and ^ 1.7 eV above. 

For all systems analyzed, we find that the essential spin 
interactions are described by the spin Hamiltonian 



(2) 



where F*'' are the isotropic magnetic couplings and a 
zero-field splitting parameters (ZFS), which is negative if 
z is an easy axis; the z direction is defined as the axis per- 
pendicular to the ring (Fig. [T|) or to the triangle (Fig. ^ . 



The coupling F* 
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- Fgg is the sum of a ferromag- 



we obtain via cLDA calculations, and a super-exchange 
term Fgg , which can be FM or antiferromagnetic (AFM). 

In Fig. [3] we show the calculated Fgg for Crg (ScP, 
= 3/2) as a function of J7 = W^' and J = J*'\ The 
figure can be understood from the analytical expression 
of Fg^ in the limit in which only density-density Coulomb 
interactions and leading order terms are retained. 
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Eq. ^ shows the competition between the first two 
terms, which yield a positive, i.e., AFM contribution and 
the FM third term, arising from excitations to empty 
states. For realistic parameters, Fgg is small and AFM. 

The ZFS term in ^ originates from the combined 
action of crystal-field and spin-orbit interactions. In the 
case of Crg it is given by 
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Here \S,M) are many-electron states in the S ~ 3/2 
ground multiplet with energy Es, while |m) are all the 
excited multiplets connected to | S, M) by the spin-orbit 
interaction, and have energy £",„. Our calculations yield 
the full ZFS tensor, and thus we can identify the easy 
magnetization axis, which in general is site-dependent. 
Remarkably, we find that the molecular global z axis (see 
above) is a nearly-easy axis for all sites. 

In the next paragraphs we discuss the results of our 
calculations (Tab. HI)) , in comparison with experiments. 
Let us start from the AFR Crg. Fig. |3] shows that, for 
realistic U and J values, the super-exchange coupling 



-'^SE^^ is AFM and of the order of few mcV. By using the 



netic (FM) screened Coulomb exchange term, F^g , which 



value of U and J obtained in cLDA (Tab. [IT| wc obtain 
Fg^"^^ — 1-99 mcV. We find that the ferromagnetic direct 
Coulomb exchange is F^^^ = —0.34 meV. Hence, the to- 
tal Heisenberg exchange constant is F''*"*"^ = 1.65 meV, 
AFM and in excellent agreement with experiments [loj . 
In addition, we find that the next-nearest neighbors ex- 
change interaction is tiny ( F*''+2 ~ io~^ \]^\^ 
explains why all experimental data can be interpreted on 
the basis of a nearest-neighbor spin Hamiltonian. Be- 
side the dominant isotropic exchange coupling, we also 
find a sizable single-ion ZFS term. Our calculations 
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Org Fe4 CrrNi 

Fgi^ 1.99 3.25 2.10 

r^;? 1.99 0.15 1.99 





Th Exp na 


Th ExpB 


Th Expia 


pi, 2 


1.65 


1.46 


2.45 2.05 


1.75 1.70 


p2,3 


1.65 


1.46 


-0.08 -0.09 


1.65 1.46 




-0.06 


-0.03 


-0.03 


-0.48 -0.35 




-0.06 


-0.03 


-0.06 


-0.06 -0.03 



TABLE II: Top: Calculated super-exchange exchange cou- 
plings. Bottom: Calculated total magnetic couphngs [3l[ 
and zero- field splitting (Th) versus experiments (Exp). Sites 
i = 1, 2, 3 are defined in Fig. [T] and Fig. [1 In CrrNi the Ni 
ion is on site 1 of Fig. [T] 

yield a significant easy-axis anisotropy in the z direction 
(D' < 0); non-axial terms are an order of magnitude 
smaller than in line with experiments. The calcu- 
lated is twice the value extracted from inelastic neu- 
tron scattering data, a remarkably good agreement given 
the small value of D*. 

Next we consider the SMM Fe4 (Fig. [J). This molecule 
has D3 symmetry; three Fe'^"'' ions (3(i^, S = 5/2) are lo- 
cated at the vertices of an equilateral triangle, and the 
fourth is at its center We find an AFM isotropic 
magnetic coupling between the central and external ions 
(r^'^) and a small FM interaction between the external 
ions (r^^'^), in excellent agreement with the values de- 
termined from experiments 0. We find that the super- 
exchange term is small for the external ions, and thus 
the FM Coulomb exchange dominates. Finally, we cal- 
culate the ZFS tensor and find again a nearly easy-axis 
anisotropy along z. Our findings for differ of less than 
a factor 1.5 from the experimental results [s^ . 

As last case we consider CrrNi, an heteronuclear AFR 
that can be obtained from Crg by replacing a Cr'^^ ion 
with a Ni2+ (3d®, 5" = 1). This system is theoretically 
the most challenging, because two different types of ion 
(Cri'^"'" and Ni'^"'") are present. Again, we reproduce well 
all experimental results. We find that the total Cr-Ni 
isotropic coupling is AFM, F^-^ = 1.75 meV, while the 
Cr-Cr coupling F^'^ is slightly smaller and close to the 
value obtained for Crg. The ZFS parameters obtained 
for Ni^+ ion [3^ are again negative in sign (easy-axis 
anisotropy) and much larger than those of Cr^"*", in agree- 
ment with neutron spectroscopy results [l^l . 

In conclusions, we present an ab-initio approach to cal- 
culate the terms of the spin Hamiltonians for molecular 
nanomagnets. It is based on the construction of many- 
body Hubbard-like models, using Foster-Boys orbitals as 
a one-electron basis [3^ . We show that this scheme works 
remarkably well for MNMs. For all systems considered, 
our results are closer to the experimental finding than 
those obtained by total-energy spin-configurations calcu- 



lations based on the B3LYP functional [l^. Differently 
than spin-configurations based approaches, our method 
allows us to determine the spin models without a priori 
assumptions on the form and the range of the Hamil- 
tonian; furthermore, since it yields the parameters of 
the Hubbard model, it works also when charge fluctu- 
ations are sizable and the spin is not well defined, like 
for molecules with metal- metal bonds [33| , or when elec- 
trons are less localized, such as in 4c? and 5d systems, 
and can be used to calculate spectra. This scheme is 
flexible, its complexity does not increase when the sym- 
metry is low, and it docs not rely on the B3LYP or 
LDA-|-[/ approximation to correlation effects. It could 
become essential for modeling MNMs whose spin Hamil- 
tonian contain many anisotropic terms, in particular if 
the principal-axis directions and relative magnitude can- 
not be inferred simply by inspecting the molecular struc- 
ture, as often the case for Co or /-electron systems. Fi- 
nally, the many-body models for small MNMs can be still 
exactly solvable, allowing to test approximations often 
adopted but impossible to test in bulk correlated sys- 
tems. Thus, we believe that our approach could become 
the method of choice for exploring fundamental issues 
and testing approximations, and for identifying and de- 
signing new molecules for quantum devices. 

Calculations were done on the Jiilich supercomputer 
Juropa, grant number JIFF46. E.P. acknowledges finan- 
cial support from the Deutsche Forschungsgemeinschaft 
through research unit FORI 346. 
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